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quark matter. The stability of strange quark matter with different strangeness 
fraction is studied. The interaction between quarks and vector mesons desta- 
bilizes the strange quark matter. If the strength of the vector coupling is the 
same as in hadronic matter, strangelets can not be formed. For the case of 
[3 equilibrium, there is no strange quark matter which can be stable against 
hadron emission even without vector meson interactions. 
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^ I. INTRODUCTION 

Strange quark matter has attracted a lot of interest since Witten suggested that it 
could be absolutely stable even at zero temperature and pressure [p]]. The investigation of 
such a possibility is relevant not only for high energy physics, but also for astrophysics. 
For example, the core of a neutron star may be composed of quark matter. The possible 
existence of strange stars which are made entirely of deconfined u, d and s quarks is one 
of the most intriguing aspects of modern astrophysics. There have been some reports of 
events with A ~ 350-500 and Z ~ 10-20 in cosmic ray experiments ||- J3J, the so-called 
exotic cosmic ray events. Also, recent studies have shown that X-ray burst sources are 
likely strange star candidates ||- M. It is also interesting to produce strange quark matter 
(strangelets) in the laboratory because they could serve as a signature of the formation of the 
quark-gluon-plasma which is a direct demonstration of QCD ||- pT| . Many ultrarelativistic 
heavy-ion collision experiments at Brookhaven and CERN |12 are proposed to search for 



X 



(meta) stable lumps of such kind of strangelets. Recently, Ardouin et al. |13|] presented a 
novel method which can be applied to characterize the possible existence of a strange quark 
matter distillation process in heavy-ion collisions. Up to now, there is no experiment which 
confirms the existence of strangelets. For example, the E864 collaboration found that there 



is no evidence for strangelet production in 11.5GeV/c per nucleon Au+Pb collisions JT4 . 

Besides the experimental efforts, there are also a lot of theoretical investigations of the 
stability of strange quark matter. The earliest discussions are based on the MIT bag model 
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|15|| which assumes that quarks are confined by a phenomenological bag. Within the bag 



quarks are asymptotically free. Calculations (jlfj within this model indicate that there is 
a range of parameters in which strange quark matter is absolutely stable, i.e. the energy 
per baryon is less than 930 MeV. The stability of strange quark matter with finite volume 
(strangelets) was also discussed in the MIT bag model. Berger and Jaffe |T7| discussed the 
surface correction for the strangelets, where they found that the surface tension destabilizes 
strangelets. The curvature contribution was considered by Madsen [P| which is dominant 
for strangelets with small baryon numbers. Though the bag model is simple, it is an in- 
complete description of confinement. Results from lattice calculations |TU| show that quark 



matter does not become asymptotically free and some hadronic degrees of freedom remain 
within the quark matter immediately after the phase transition. Fowler, Raha and Weiner 
PD]| suggested another description of the confinement mechanism via the introduction of 



a density-dependent quark mass. This quark mass-density-dependent (QMDD) model was 
first employed to study the properties of ordinary quark matter [^TJ and then applied to the 
investigation of strange quark matter pl| |- p5| |. As was pointed in our recent paper [|26 



their thermodynamic treatment was not correct. We reconsidered strange quark matter in 
the self-consistent quark mass density dependent model and found a region of parameters 
in which the strange quark matter is absolutely stable. 

In the QMDD model, the concept of a density dependent quark mass has no dynamical 
origin. In recent years, some approaches for strange quark matter based on dynamical models 
were developed. Alberico et al. |27[ utilized the color dielectric model to calculate the energy 
per baryon of strange quark matter. They found that while the double minimum version of 
the color dielectric model allowed the existence of strangelets, the single minimum version of 
this model excluded the possibility. Stability of strange quark matter was also investigated 
using the effective 4-quark interactions pg| , the SU(3) Nambu-Jona-Lasinio (NJL) model 
with and without 4-quark vector type interactions p9| , |30l . In studying hadronic matter, we 
proposed a chiral SU(3) quark mean field model. This chiral quark model was applied to 
investigate the properties of strange hadronic matter and multi-strange hadronic systems 



3l| , |32|| . In this paper, we want to use this model to discuss the stability of strange quark 



matter. The difference between quark and hadronic matter is that in quark matter, the u, 
d, s quarks are deconfined and not combined into baryons by the confining potential. 

The paper is organized as follows. In section II, we introduce the basic model features. 
We apply the model to investigate strange quark matter in section III. The numerical cal- 
culations are discussed in section IV. Finally, main conclusions are drawn in section V. 



II. THE MODEL 

Our considerations are based on the chiral SU(3) quark mean field model (for details 



see Refs. [^T|,^|). For completeness, we introduce the main concepts of the model in this 
section. In the chiral limit, the quark field q can be split into left and right-handed parts qi 
and qR. q = qi + qn- Under SU(3)l x SU(3)r they transform as 

q' L = L( 1l, q' R = Rqn- (1) 
The spin-0 mesons are written in the compact form 
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M(M+) = E ± ffl = —= (cr a ± ivr a ) A a , 

V 2 a=0 



(2) 



where a a and 7r a are the nonets of scalar and pseudoscalar mesons, respectively, X a (a = 
1, 8) are the Gell-Mann matrices, and A = ^/| /. Plus and minus sign correspond to M 
and M + . Under chiral SU(3) transformations, M and M + transform as M — > M' = LMR + 
and M + — > M + = RM + L + . As for the spin-0 mesons, the spin-1 mesons are set up in a 
similar way as 



l,{r,)= l -{V,±A,) 



2^=0 



(3) 



Rr fl R + . These matrices can be written in a form 



They transform as / M 

where the physical states are explicit. For the scalar and vector nonets, the expressions are 
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Pseudoscalar and pseudovector nonet mesons can be written in the same way. 

The total effective Lagrangian for the description of strange quark matter is given by: 
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It contains the free part for massless quarks C q o = q i^f^d^ q, the quark-meson field interaction 
term 

£qM = 9s (<iLMq R + q R M + qi)j - g v (q L Y^QL + <?i?7%^) , (7) 
the chiral-invariant scalar meson and vector meson £yy self-interaction terms in the 



mean field approximation [pT| , |33 



Cm = -\ hx 2 2 + C 2 ) + h (a 2 + c 2 ) 2 + h + c 4 ) + h X a 2 ( 



1 x 4 $ 
-hx 4 - ~X 4 ln— + ^X 4 ln- , 
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£vv = ^^ (mluj 2 + m 2 p p 2 + m 2 2 ) + g 4 {to 4 + 6cu 2 p 2 + p 4 + 



(8) 
(9) 



where 6 = 6/33; a , Co an d Xo are the vacuum values of the mean fields a, ( and x an d 
the three terms £ x sb, ^Am, and Ch which explicitly break the chiral symmetry. Chiral 
symmetry requires the following basic relations for the quark-meson coupling constants: 
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(10) 
(11) 



Note, the values of a , ( and xo are determined later from Eqs. (P0|) - (p2|) . Particularly, the 
parameters <To and Co are expressed through the pion (F w = 93 MeV) and the kaon (Fx = 
115 MeV) leptonic decay constants as: 
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The Lagrangian C x sb generates the nonvanishing masses of pseudoscalar mesons 
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leading to a nonvanishing divergence of the axial currents which satisfy the PCAC relations 
for 7T and K mesons. Scalar mesons obtain the masses by spontaneous breaking of the 
chiral symmetry in the Lagrangian (§). The masses of u, d and s quarks are generated by 
the vacuum expectation values of the two scalar mesons a and C- To obtain the correct 
constituent mass of the strange quark, an additional mass term should be added: 



—Am s qSq 



(14) 



where S = | (i — Asv^) = diag(0, 0, 1) is the strangeness quark matrix. Finally, the quark 
masses are given by 



m u = m d = -—^vo, rn s 



-g s Co + Am s , 
29 MeV are determined from m„ 



(15) 



The parameters g s = 4.76 and Am s = 29 MeV are determined from m q = 313 MeV and 
m s = 490 MeV. In order to obtain reasonable hyperon potentials in hadronic matter we 
include the additional coupling between strange quarks and scalar mesons a and C [|3l 
This term is expressed as 



£h = (hi a + h 2 C) ss . 



(16) 



III. APPLICATION TO STRANGE QUARK MATTER 

Now we apply the model to investigate strange quark matter. We begin with the ther- 
modynamical potential because all other quantities such as energy per volume and pressure 
can be obtained from it. The thermodynamical potential is defined as 

q= 53 ~ 2 ( ^p 7r r Sk { ln i 1 + e~^-^ kBT ) + In (l + e -re(*)+M/fc*r)| _ Cm , (17) 



=q,e 
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where E*(k) = ^Jm* 2 + k 2 , 7 T is 3 for quarks and 1 for electrons and Cm is the meson 
interaction including the scalar meson self-interaction £es, the vector meson self-interaction 
Cyy and the explicit chiral symmetry breaking term C x sb- I n the MIT bag and the QMDD 
model, Cm is replaced by the effective bag constant. At zero temperature, fl can be expressed 
as 

1 



n = - £ 
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+ 3m* 4 ln 
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(18) 



where /x e is the chemical potential of the electron and the quantity z/j (i—u, d, s) is related 
to the usual chemical potential fa by Vi — fa — g^u — g^<p. The effective quark mass is given 
by m* = —g l a cr — g^( + m i0 . The total baryon density is defined as 



PB = + Pd + Pa)- 



(19) 



With the thermodynamical potential, the energy per volume e and pressure p of the system 
can be derived as e = Q + J2i= u ,d, s ,e PiPi an d P — ~ ^- 

The mean field equations for the meson 0j are obtained with ^ = 0. For the scalar 
mesons a, ( and x, the equations are expressed as 
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The equations for the vector mesons can be obtained in the same way as 
X 2 
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The scalar and vector densities can be written as 

(2tt)3 7 ^k 2 + mf 

<^>=7^t^ (27) 



with /c^ = Juf — m* 2 . 

IV. NUMERICAL RESULTS 

Now we investigate the implications of the previous formalism to strange quark matter. 
Compared to earlier applications to hadronic matter, here we need not to introduce the 
confining potential which combines quarks into baryons. In quark matter, the u, d and s 
quarks are deconfined. They only interact by scalar and vector mesons. The self-interactions 
between mesons are the same as in hadronic matter. All the parameters in this model 
are determined in our previous papers. They are listed in Table I. We assume that these 
parameters do not change when the model is applied to quark matter. In fact, the parameters 
which describe the interactions between fields 'should' be universal. The medium effects are 
included by the treatment (relativistic mean field theory) itself and are not included in the 
original Lagrangian. In our discussion, we first do not include the additional coupling Ch 
and discuss the effect of this term on strange quark matter later. 

If strange quark matter is stable and can survive for a long time, equilibrium with respect 
to the weak processes: 

s — > u 4~ e~ 4~ u e , d — > u + e~ + u e , (28) 

may be achieved. If we neglect the chemical potential of the neutrino, the chemical potentials 
of quarks and electron have the following relation: 

= Us = + (29) 

The values for p u and fid are determined by the baryon density and the total charge Q. As 
usually, we assume the charge Q to be zero. The equations for the mesons (^0|)-(^5|) can be 
solved simultaneously. The energy per baryon with different strength of vector coupling and 
for different parameter sets are shown in Fig.l. When vector interactions are not considered, 
there is a local minimum of energy per baryon with parameter set A. The corresponding 
density is about 0.13/m~ 3 . As shown later, at this density no strange quarks appear. When 
the vector interactions are included, even when the strength is half of that in hadronic 
matter, the local minimum for the energy per baryon disappears. The energy per baryon 
E/A increases monotonously with the increasing baryon density. 

In Fig. 2, we show the fractions = of u, d, and s quark versus baryon density with 
and without vector interactions. The fraction of u quarks almost does not change with the 
density. There exists a relationship p d + p s ~ 2p u . Therefore, although we included the 
electron in our calculations, the fraction of electrons is very small. The chemical potential 
of the electron p e is not zero which results in a larger fraction of d than u quarks. At low 
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density no strange quarks are present. The fraction of d quarks is about two times that 
of u quarks. Without vector interactions, when the density is larger than about 0.41/m~ 3 , 
strange quarks appear. Compared to Fig.l, this density is larger than the density where 
the local minimum appears. This result is close to that of Ref. J3(| where the SU(3) NJL 
model is used. Therefore, no stable strange quark matter can exist at zero pressure in 
contradiction to the original suggestion by Witten [p]] The so-called strange star cannot be 
composed entirely by the deconfined u, d and s quarks. The stable strange quark matter 
can only exist in the core of these objects where the pressure is high enough to force the 
transition from hadronic matter to quark matter to occur. 

If strange quark matter is metastable, then it may be produced in heavy ion collisions. 
In this case, the /3-equilibrium may not be achieved. In our calculations, we assume that 
/J-u = Hd = /V The values of \i q and ji s are determined by the total baryon density and the 
strangeness fraction f s (f s = 3r s ). In Fig. 3, we plot the effective masses of nonstrange and 
strange quarks versus the baryon density for different strangeness fraction with parameter 
set A. Both u (d) and s quark masses decrease with the increasing density. In nonstrange 
quark matter, the mass of the u (d) quark decreases more quickly than that of the s quark. 
When f s increases, the mass of the s quark decreases at fixed density. At some high density, 
the strange quark mass is even lower than the mass of nonstrange quarks. Compared to 
the QMDD model, here the effective quark masses are obtained dynamically. The quark 
masses are not only density dependent but also strangeness fraction dependent which are 
not present in the QMDD model. 

In Figs. 4-6 we plot the energy per baryon versus density for different values of f s , which 
corresponds to different vector coupling constants. The solid and dashed lines are for pa- 
rameter sets A and B, respectively. First we do not include the interactions between quarks 
and vector mesons. This is close to the QMDD model or SU(3) NJL model with only scalar 
type 4-quark interactions. In Fig.4, for parameter set A, there exists a local minimum of 
energy per baryon for any f s . The baryon density at the minimum of energy per baryon first 
increases and then decreases when f s increases. For nonstrange quark matter, though there 
is a local minimum, the system has no positive binding energy compared to the vacuum 
mass of nonstrange quarks. For strange quark matter, the energy per baryon is lower than 
the masses of hyperons with the same strangeness number. When f s — 1, the binding energy 
is about 45 MeV compared to the vacuum constituent quark mass. The maximum binding 
energy is about 60 MeV and the corresponding f s is about 2.0. Metastable strange quark 
matter is therefore favored to have a large strangeness fraction (high negative charge). For 
the parameter set B, the system has a local minimum of energy per baryon only for some 
range of f s , that is 1 < f s < 3. The maximum binding energy is about 5 MeV, when f s is 
around 2.0. 

When the interactions between quarks and vector mesons are included, the system is 
destabilized. In Fig. 5, we plot the energy per baryon versus baryon density with the vector 
coupling constant half of that in hadronic matter. Compared to Fig.4, the energy per baryon 
becomes higher. For parameter set A, there exists a local minimum for < f s < 2. The 
maximum binding energy is only about 10 MeV with f s ~ 1.0. The corresponding density is 
also much lower compared to Fig.4. For parameter set B, there is no local minimum for any 
strangeness fraction. If the vector coupling is the same as in hadronic matter, as is shown 
in Fig. 6, the energy per baryon will increase monotonously with the baryon density for both 
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parameter sets A and B. The vector meson couplings are also discussed in Ref. | 29| . Though 
these two models are quite different, the results are comparable to each other. 

If vector interactions are included, even when the strength is only half of that in hadronic 
matter, the binding energy of metastable strange quark matter becomes small or negative. 
At the minimum of the energy per baryon, the pressure p of the system is zero. This 
kind of objects with finite volume are called strangelets. If the vector interactions are fully 
considered, the metastable strangelets can not be formed at zero pressure. 

Up to now, we did not consider the additional coupling between strange quark and scalar 
mesons Ch, which is important to obtain reasonable hyperon potentials in hadronic matter 
pl| . When the additional term is included, the effective quark masses versus baryon density 
for parameter set A are given in Fig. 7. The result for the strange quark mass evidently 
change, especially for large strangeness fraction, when compared to Fig. 3. In nonstrange 
quark matter, the effective mass of s quarks almost stays constant. When the strangeness 
fraction is high, the strange quark has a lower mass when compared to Fig.3, where Ch is 
not included 

The effective quark masses will affect the energy of the system. In Fig. 8 we plot the energy 
per baryon versus density with the vector coupling constant half of that in hadronic matter. 
The solid and dashed lines correspond to the cases without and with C h , respectively. For 
small strangeness fraction f s , the results of these two cases are close. For large strangeness 
fraction, the additional term Ch will produce larger binding energy. If the vector coupling 
strength is as big as in hadronic matter, there is no local minimum for both cases, with 
or without Ch- Therefore, the inclusion of the additional term Ch does not affect the main 
results of Fig.l. This is because when the /3-equilibrium is achieved, the strangeness fraction 
of the system is smaller than 1. The additional term Ch only gives sizable contributions for 
systems with a large strangeness fraction. 



V. CONCLUSIONS 

We investigate strange quark matter in a chiral SU(3) quark mean field model. The 
effective quark masses are obtained dynamically by the quark meson interactions and they 
are both density and strangeness fraction dependent. The stability of strange quark matter 
is studied for different values of the vector coupling constant and for different parameter 
sets. The effect of the additional term Ch is also discussed. 

If the strange quark is stable and can survive a long time, the /3-equilibrium can be 
achieved. The strangeness fraction f s is smaller than 1. In the chiral SU(3) quark model, 
at the density where the system has a local minimum for the energy per baryon, no strange 
quarks appear. Even without the vector meson coupling, this nonstrange quark matter 
has a negative binding energy and cannot be bound. Therefore, opposed to the suggestion 
by Witten, stable quark matter at zero pressure cannot exist even without vector meson 
interactions. 

If quark matter can be produced in heavy ion collisions, the /3-equilibrium may not be 
achieved. This metastable strange quark matter can have a high strangeness fraction (high 
negative charge). When we do not take the vector meson interaction into account, as in the 
QMDD model or SU(3) NJL model with only scalar type 4-quark interaction, the maximum 
binding energy is about 5 MeV - 60 MeV and the corresponding strangeness fraction is 
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about 1-2. If we assume that the strength of vector meson coupling is half of that in 
hadronic matter (in this case, the strength of vector coupling is comparable to the one of 
Ref. [29]), the maximum binding energy decreases. Inclusion of the additional term Ch will 
forces the system to have a larger strangeness fraction and binding energy. When the vector 
coupling constant is the same as for hadronic matter, no metastable strange quark matter 
(strangelets) can be formed in heavy ion collisions. 

However, strange quark matter can still exist in the core of a neutron star. The phase 
transition from hadronic matter to quark matter can occur at high density and pressure 
which is caused by gravity. Because both hadronic and quark matter can be described in 
the SU(3) quark mean field model, it is of interest to investigate this phase transition in this 
model. This will be studied in the future. 
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TABLES 

TABLE I. Parameters of sets A and B of the model. 



set 


k 


h 


k 2 


h 


&4 


9s 


9v 


9a 




h 2 


A 


4.94 


2.12 


-10.16 


-5.38 


-0.06 


4.76 


10.92 


37.5 


-2.20 


3.24 


B 


3.83 


2.64 


-10.16 


-3.40 


-0.18 


4.76 


10.13 


0.0 


-2.03 


2.55 



Figure captions 

Fig.l: Energy per baryon versus baryon density for different vector coupling constants. 
The solid and dashed lines are for parameter sets A and B. Results are shown for the case 
of /^-equilibrium. 

Fig. 2: Fractions of u, d and s quarks in strange matter versus baryon density with and 
without vector interactions, respectively. The solid and dashed lines are for parameter sets 
A and B. Results are shown for the case of /5-equilibrium. 

Fig. 3: The effective nonstrange and strange quark mass versus baryon density for dif- 
ferent strangeness fraction. Parameters are taken from set A. 

Fig. 4: The energy per baryon versus baryon density for different strangeness fraction for 
parameter sets A and B. The interaction between quarks and vector mesons is not considered. 

Fig. 5: Same as in Fig. 4, but here the quark- vector meson interaction is included. The 
value of the vector coupling constant is half of that used in hadronic matter. 

Fig. 6: Same as Fig. 5, but the vector coupling constant is of the strength as in hadronic 
matter. 

Fig. 7: Same as Fig. 3, but the the additional term C h is included in the calculation. 

Fig. 8: The energy per baryon versus baryon density calculated with parameter set A. 
The vector coupling constant is half of that in hadronic matter. Solid and dashed lines 
correspond to the case without and with the additional term Ch, respectively. 
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